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Summary. — We provide measurements of the neutral hydrogen fraction xm at 
z ^ Q,hj comparing semi-analytical models of the Lya forest with observations of 
high-2 quasars and Gamma Ray Bursts absorption spectra. 

We analyze the transmitted flux in a sample of 17 QSOs spectra at 5.74 < Zem < 6.42 
studying separately the narrow transmission windows (peaks) and the wide dark 
portions (gaps) in the observed absorption spectra. By comparing the statistics 
of these spectral features with our models, we conclude that xm evolves smoothly 



from 10 



at z 



5.6, with a robust upper limit xwi < 0.36 at 



z — 6.3. We show the results of the first-ever detected transverse proximity effect in 
the HI Lya forest, produced by the HIT region of the faint quasar RD J1148-I-5253 
at z = 5.70 intervening along the LOS of SDSS J1148-h5251 at z = 6.42. 
Moreover, we propose a novel method to study cosmic reionization using absorption 
line spectra of high-redshift GRBs afterglows. We show that the time evolution and 
the statistics of gaps in the observed spectra represent exquisite tools to discriminate 
among different reionization models. By applying our methods to GRB 050904 
detected at z = 6.29, we show that the observation of this burst provides 
strong indications of a highly ionized intergalactic medium at z ~ 6, with an esti- 
mated mean neutral hydrogen fraction xui ~ 6.4±0.3 x 10~^ along that line of sight. 
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1. Introduction 

Although observations of cosmic epochs closer to the present have indisputably shown 
that the InterGalactic Medium (IGM) is in an ionized state, it is yet unclear when the 
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phase transition from the neutral state to the ionized one started. Thus, the redshift of 
reionization, Zrei, is still very uncertain. 

In the last few years, our knowledge of the reionization process has been enormously 
increased mainly owing to the observation of 2 ~ 6 QSOs by the SDSS survey [2] and 
CMB data [9] . Long gamma ray bursts (GRB) may constitute a complementary way to 
study the reionization process possibly probing even larger redshifts, the current recold 
holder being GRB 050904 at ^ = 6.3 [10] [6]. 

We provide measurements of the neutral hydrogen fraction xm at epochs approaching 
reionization, by comparing semi-analytical models of the Lya forest with observations of 
the highest-z QSOs and GRBs absorption spectra. 

2. — Lya forest simulation 

The ultraviolet radiation emitted by a QSO/GRB can suffer resonant Lya scattering 

as it propagates through the intergalactic neutral hydrogen. In this process, photons 
are removed from the line of sight (LOS) resulting in an attenuation of the source flux, 
the so-called Gunn-Peterson (GP) effect. To simulate the GP optical depth (top) dis- 
tribution we use the method described by Gallcrani et al. (2006) [3] and further revised 
in Gallerani et al. (2007) [4], whose main features are summarized as follows. Mildly 
non-linear density fluctuations giving raise to spectral absorption features in the inter- 
galactic medium (IGM) can be described by a Log-Normal distribution. For a given 
IGM equation of state, the mean HI fraction, xhi, can be computed from photoioniza- 
tion equilibrium as a function of the baryonic ovcrdensity, A = p/p, and photoionization 
rate, F, due to the ultraviolet background radiation field. These quantities must be de- 
termined from a combination of theory and observations; here we follow the approach of 
Choudhury & Fcrrara (2006) [1]. This model contains two free parameters: (i) the star- 
formation efficiency /*, and (ii) the escape fraction fesc of ionizing photons from galaxies. 
These are calibrated by a maximum-likelihood procedure to a broad observational data 
set. Currently, the available data can be explained by two different reionization histo- 
ries, corresponding to different choices of the free parameters: (i) an Early Reionization 
Model (ERM) (/* 0.1; fesc = 0.07), characterized by a highly ionized IGM at 2; > 6, 
and (ii) a Late Reionization Model (LRM) (/* = 0.08; fesc = 0.04), in which reionization 
occurs at z « 6. Both ERM and LRM by construction provide an excellent fit to the 
mean neutral hydrogen fraction evolution experimentally deduced from the GP test. 

3. — QSOs absorption spectra 

We first test the predictions of our models by applying various statistical analysis 
to the simulated spectra and comparing our results with observations. Specifically, we 
use the following control statistics: (i) Mean Transmitted Flux evolution in the redshift 
range 2 — 6; (ii) Probability Distribution Function (PDF) of the transmitted flux at the 
mean redshifts z = 5.5,5.7,6.0; (iii) Gap Width (GW) distribution in 3.5 < 2: < 5.5. For 
what concerns the GW statistics we define gaps as contiguous regions of the spectrum 
having a tgp > 2.5 over rest-frame wavelength (Xrf) intervals > 1 A. As both ERM 
and LRM successfully match the observational data at 2; < 6 for the control statistics 
considered, we proceed the comparison with more advanced statistical tools. 

3'1. Largest Gap Width analysis. - The Largest Gap Width (LGW) distribution quan- 
tifies the fraction of LOS which are characterized by the largest gap of a given width. 
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Fig. 1. - Largest Gap Width distribution for the LR and the HR cases (left and right, respec- 
tively). Filled circles represent the result of the analysis of the 17 QSOs observed spectra. Solid 
red (dotted blue) lines show the results obtained by the semi- analytical modeling implemented 
for the ERM (LRM). Vertical error bars measure poissonian noise, horizontal errors define the 
bin for the gap widths. 



We use observational data including 17 QSOs obtained by Fan at al. (2006). We divide 
the observed spectra into two redshift-selected sub-samples: the "Low-Redshift" (LR) 
sample (5.7 < z^^. < 6), and the "High-Redshift" (HR) one (6 < Zem < 6.4). In order 
to measure the evolution of xhi with redshift, we apply the LGW both to simulated and 
observed spectra. From the comparison shown in Fig. 1 it results that both the LGW 
distributions predicted by the ERM and LRM provide a good fit to observational data. 
We exploit the agreement between the simulated and observed LGW distributions to 
derive an estimate of xm- We find log;^Q xm = — 4.41q;9q at z ~ 5.3(^). By applying the 
same method to the HR sample we constrain the neutral hydrogen fraction to be within 
logioa^Hi = -4.2^?:^^ at z«5.6. 

Although the predicted LGW distributions are quite similar for the two models consid- 
ered, yet some differences can be pointed out. Indeed in the HR case we find that a 
neutral hydrogen fraction at z « 6 higher than that one predicted by the LRM would 
imply an even worst agreement with observations, since a more abundant HI would pro- 
duce a lower (higher) fraction of LOS characterized by the largest gap smaller (higher) 
than 40 A with respect to observations. Thus, this study suggests xm < 0.36 at z — 6.32 
(obtained from the maximum value for xhi found in the LRM at this epoch). 

3'2. Transverse proximity effect. - Transmissivity windows in absorption spectra can 
be produced by ionizing sources whose bubbles intersect the lines of sight to the target 
object. Mahabal et al. (2005) [8] have discovered a faint quasar (RD J1148-I-5253, 
hereafter QSOl) at z = 5.70 in the field of the highest redshift quasar currently known 
(SDSS J1148-I-5251, hereafter QS02) at z = 6.42. In this Section we study the QS02 
transmitted flux, in order to analyze the proximity effect of QSOl on the QS02 spectrum. 

In Fig. 2 (left panel) we compare the observed transmitted flux in the spectrum of 



(^) The Xm value quoted is the mean between the estimates predicted by the ERM and the 
LRM. Moreover, we consider the most conservative case in which the errors for the measurement 
of the neutral hydrogen fraction are provided by the minimum xm value found in the ERM and 
the maximum one in the LRM. 
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Fig. 2. - Left Leftmost panels: Observed transmitted flux (black spectra) in the spectrum of 
QSO SDSS J1148+5251 (QS02, z^m = 6.42). The solid black line shows the redshift path 
(Azprox) in which the bubble produced by QSO RD J1148+5252 (QSOl, Zem = 5.65) inter- 
sects the LOS to QS02. Top panels (ai), with i=l,3: Simulated fluxes (cyan spectra) along 
3 different random LOS (cases "without bubble"). Bottom panels (bi), with i—1,3: Simulated 
fluxes (magenta spectra) along the same LOSs shown in the top panels, taking into account 
the contribution from QSOl to the total ionizing flux (cases "with bubble"). Right Evolution 
of the optical depth r as a function of the distance R from QSOl. Filled circles denote the 
observed mean value for r, while error bars represent the maximum and the minimum observed 
T at a given distance from the foreground QSO. Solid (dotted) magenta lines are the mean 
(maximum/minimum) values from 500 simulated LOS, computed adopting the case "with bub- 
ble" . The dashed cyan horizontal line shows the mean optical depth predicted by the ERM in 
correspondence of the emission redshift of the foreground QSO. The dotted cyan horizontal lines 
denote the maximum/minimum optical depth at the same redshift. 



QS02 with the simulated fluxes along 3 different LOS with (bottom row) or without (top) 
including the contribution from QSOl to the total ionizing flux. For brevity, we refer to 
these case as "with bubble" or "without bubble" . Visual inspection of Fig. 2 shows that 
the case "with bubble" is in better agreement with observations. Such statement can be 
made more quantitative by introducing a quantity denoted Peak Spectral Density (PSD), 
i.e. the number of peaks per unit Xrf interval. For both the observed and simulated 
spectra, we compute the PSD inside and outside the bubble, finding the following results: 

(PSDOUT^PSDi^J = (0.11,0.46); {PSDOU^^ ^ PSDli'J = (0.04t« °t 0.22t«:?2). 
Both in observations and simulations, the PSD is found to be « 4.5 times larger inside 
that bubble than outside it. 

As a final test for our model, we compute the observed evolution of the optical depth 
as a function of the distance R from QSOl and compare it with the predictions of model 
"with bubble"; the result is shown in Fig. 2 (right panel). The agreement between 
observations and simulations is at 1-a confidence level for 70% of the plotted points. For 
R < 2 Mpc, the mean optical depth 1.5 < f < 3.5 is lower than the mean value expected 
at z = 5.65 (fs.es 4); it approaches fs.gs at distances larger than Rr ^ 2 Mpc. By 
taking the difference between Rt and R±, we set a lower limit on the foreground QSO 
lifetime tg > ^^^^^-^ _(- (t^ _ igsoi) ~ H Myr, where tr and tqsoi represent the cosmic 
times corresponding to the redshifts Zr — 5.68 and z^^f^^ = 5.65, respectively. 

The physical interpretation of the results reported in this Section is the following. In 
the case "with bubble", in correspondence of spectral regions where r^f*^^ > Fhi, most 
of the gaps present in the case "without bubble" disappear, making room for peaks, as 
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a consequence of the decreased opacity in the proximity of QSOl. The enhancement in 
the transmissivity decreases moving toward outside the bubble, as T^^'^^ < Fhi- This is 
the first-ever detection of the transverse proximity effect in the HI Lya forest. 

4. GRBs absorption spectra 

We have built a database of synthetic GRB afterglow emission spectra starting from 
the observed spectral energy distribution and time evolution of the most distant GRB 
detected up-to-now, i.e. GRB 050904 [10] [7]. The unabsorbed afterglow spectrum of 
GRB 050904 can be parameterized as F{v) cx v'^t^ , where we assume for a and /? the 
values founded by previous studies [10] [6]. Finally, we normalize the intrinsic GRB 
050904 optical spectrum in order to reproduce the flux of ~ 18 /xJy as measured at 1 
day from burst in the J band [6]. We simulate the observed flux Fobs = ^(i') e~'^°^ of 
GRBs absorption spectra in the rest frame spectrum between Lya (1215.67 A) and Ly/3 
(1025.72 A) at different times after the burst. 



150 


' 1 


,=6.3 

A 




,,,,,,,,,,,, 1 _ 

■_ A y : 


100 








■ f" ■ 


50 


- i 

4 






: 1 : 





" . 1 


. 1 . , . 1 . . 




" , 1 , , , 1 , , , 1 , , , " 





1 2 3 4 5 
time[d] 




1 2 3 4 5 
time[d] 



Fig. 3. - Evolution of the largest transmissivity gap found in the rest frame GRB afterglow 
spectra as a function of the observer time after the explosion for the two selected GRB redshifts 
z — 6.3 (left panel) and z = 6.7 (right). The red solid (blue dashed) line refers to Early 
(Late) Reionization Model, with the error bars showing the standard deviation associated to 
the weighted mean computed from 10 realizations of 5 LOS. The black filled circle refers to the 
largest gap measured in the spectrum of GRB 050904 afterglow. 



4'1. Largest Gap Width analysis. - The main idea we propose is to exploit the statistics 
of the transmissivity gaps imprinted by the intervening IGM neutral hydrogen on the 
otherwise smooth power-law spectrum of high-redshift GRBs. On general grounds, we 
expect that at any given redshift, but particularly above z — 6, where differences become 
more marked, the value of xhi is higher in the LRM than in the ERM. As a result, wider 
and more numerous gaps are expected if reionization completes later. Moreover, as the 
time after the burst increases, the gaps become larger. In fact, the progressive fading of 
the unabsorbed afterglow produces a corresponding attenuation of the observed flux. 
To put the above arguments on more quantitative grounds, we have derived, using the 
procedure described above, the evolution of the LGW found in synthetic afterglow spectra 
with time after explosion. The results of the calculation are shown in Fig. 3 for two 
selected GRB redshifts zgrb = 6.3, 6.7 and for the ERM and LRM cases. The differences 
caused by the two different reionization histories are striking. Since the beginning, the 
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width of the gaps is a factor « 2 times wider in the LRM (38 A vs 24 A at 2;grb = 6.3, 
74 A vs 40 A at ^;grb = 6.7) than in the ERM. We compare these predictions to the 
results of the analysis of GRB 050904 spectrum obtained 3.4 days after the burst (Kawai 
et al. 2006), whose largest dark gap is VFmax ~ 65 A in the source rest frame. This value 
refers to the dark region immediately blueward the Lya emission line (~ [8420; 8880]). 
The black filled circle in Fig. 3, corresponding to the values of T4^max and t for GRB 
050904, clearly shows that the ERM is favored by the data. In the LRM the typical 
Wmax at t = 3.4 days is as large as ~ 100 A, well above the observed value. 
We derive the mean neutral hydrogen fraction xm along the synthetic lines of sight 
characterized by a LOW in the interval 65 ± 5 A, centered on the size of the largest dark 
gap observed in the GRB 050904 spectrum. The observed LGW in the GRB 050904 
afterglow spectrum is consistent with xm = 6.4 ± 0.3 x 10~^ [5]. This result is in 
agreement with previous measurements by Totani et al. (2006), who find that xm is 
consistent with zero with upper limit Xm < 0.17 at 68% C.L. 

5. Conclusions 

We measure the neutral hydrogen fraction xhi at epochs approaching the reionization, 
by comparing semi-analytical models of the Lya forest with observations of high-2 quasars 
and Gamma Ray Bursts absorption spectra. We consider an Early Reionization Model 
(ERM), characterized by a highly ionized Universe at 2 ~ 6 and a Late Reionization 
Model (LRM) in which reionization occurs at 2; ~ 6. 

By comparing statistical analysis of the transmitted flux in a sample of 17 QSOs spectra 
at 5.74 < Zem < 6.42 with our models, we find that both ERM and LRM provide good 
fits to the observed LGW distribution, favoring a scenario in which xhi smoothly evolves 
from 10~'* '* at z w 5.3 to 10~^'^ at z « 5.6, with a robust upper limit .thi < 0.36 at 
z = 6.3. Discriminating among the two reionization scenarios would require a sample of 
QSO at even higher redshifts. 

We show the results of the first-ever detected transverse proximity effect in the HI Lya 
forest, produced by the HII region of the faint quasar RD J1148-I-5253 at z = 5.70 
intervening along the LOS of SDSS Jl 148+5251 at z = 6.42. 

Moreover, we show that the time evolution of gaps in GRBs absorption spectra represent 
exquisite tools to discriminate among different reionization models. By applying our 
methods to GRB 050904 detected at z = 6.29, we show that the observation of this burst 
provides strong indications of a highly ionized intergalactic medium at z ^ 6, with an 
estimated mean neutral hydrogen fraction xm = 6.4 ±0.3 x 10~^ along that line of sight. 
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